Character of ground state of an aperiodic frustrated Josephson junction array 
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We study the energy spectrum for an aperiodic Josephson junction ladder, as a function of frus- 
tration. Frustration is brought about by application of a transverse magnetic field, and aperiodicity 
is imposed by the arrangement of plaquettes with two incommensurate areas. We study the effect 
of the incommensurate plaquette areas in conjunction with that of the aperiodicity. The structure 
of the energy spectrum at deep minima is shown to be described by a model that treats the pla- 
quettes independently. The energy spectrum is a quasiperiodic function of frustration; short range 
correlations in the arrangement of plaquettes have a small effect on the energy power spectrum. 
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I. INTRODUCTION 

The electronic and transport properties of quasiperi- 
odic systems are affected by their quasicrystalline struc- 
ture. In one-dimensional (ID) systems with Fibonacci or- 
der, the energy spectrum is singularly continuous and the 
wave functions are neither localized nor extended.- Simi- 
lar exotic electronic states are predicted for the 3D icosa- 
hedral quasicrystals^ The structurally ordered icosahe- 
dral phases and their crystalline periodic approximants 
show the same transport properties, implying that the 
local configurational order (common in both) plays the 
major role^ At the same time, study of the Fibonacci 
semiconductor superlatticcs indicates that at low enough 
temperatures, the longitudinal magnetoresistance of the 
quasiperiodic structure can be distinguished from that of 
its low order periodic approximant^ 

In two dimensions, two types of study have been illu- 
minating. One has been the study of the electronic struc- 
ture, the energy spectrum and wave functions, which 
again has shown results similar as above, particularly in 
case of the two-dimensional (2D) Penrose latticed The 
other kind is inspired by the Little-Parks effect^ reveal- 
ing the interplay of flux quantization and free energy and 
manifest in the variation of critical temperature of the 2D 
lattice as a function of the applied magnetic field per- 
pendicular to its p\a.neZ&- Here, we consider the ground 
state of an aperiodic array of Josephson junctions, mainly 
drawing from the second type of study. 

Josephson arrays give a useful context in which to 
study frustration. These are composed of weakly coupled 
superconductors, with an interaction between neighbor- 
ing sites, 
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integer n is chosen so that — it < Aj 



+ 2n7r, where the 
< 7r, it is readily 



H = -Ej J2 cos & 



0, - A,. 
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where 9 is the phase of the superconducting order param- 
eter on the site. Aij describes the effect of the magnetic 
field; it is proportional to the line integral of the compo- 
nent of the vector potential along the link from i to j. It 
causes the phase of the superconducting wavefunction to 
vary from place to place. If we define the gauge invariant 



seen that the Ay cannot all be zero, because the directed 
sum of the Aij around a plaquette is the magnetic flux 
through it. Further, it is found that the directed sum of 
the Ay around a plaquette sometimes differs from the 
magnetic flux by a multiple of 2ir; we will say that these 
plaquettes contain a vortex. The ground state is speci- 
fied by a configuration of vortices^ Thus these systems 
exhibit frustration^, meaning that in the ground state, 
the individual terms of Eq. (jl.ll) are not simultaneously 
at minimum. 

For a structure such as the Penrose lattice that is com- 
posed of two types of plaquettes with irrational area ra- 
tios, the plaquettes can never achieve their energy min- 
ima simultaneously as the magnetic field is varied, a phe- 
nomenon known as geometrical frustration}-- The degree 
of frustration can be characterized by a parameter /, 
which is the ratio of the external flux in a plaquette to 
the quantum of flux. Only when / is an integer for all 
plaquettes can all junctions be simultaneously at mini- 
mum energy. 

When the array is a periodic ladder and / is ratio- 
nal, the ground state is a superlattice of vortices. Vallat 
and Becki^ showed that for the cosine interaction the en- 
ergy E(f) is a uniformly continuous function of /. The 
ground state energy E(f) and the transition tempera- 
ture T c (f) are both periodic functions of /, with distinct 
local minima where / is a simple rational, showing that 
commensuration is important for this caseJ^ 

Behrooz et alZ considered the problem of flux quan- 
tization on a quasiperiodic network of thin aluminum 
wires. The network consisted of parallel evenly spaced 
wires in one direction, but with unevenly spaced paral- 
lel wires in the other. The latter spacing was generated 
by a substitution rule, so that it had inflation symme- 
try characterized by the irrational number a = 1 + v2. 
The ratio of the areas of the two kinds of plaquettes in 
the quasicrystalline network was also a. They found that 
T c (/) was no longer a periodic function of the frustration, 
but it did show a series of high peaks which in magnitude 
were nearly the same as that of the unfrustrated network. 

We consider aperiodic ladders consisting of two types 
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of plaquettes with incommensurate areas, as in Fig. 1, 
and study the ground state energy spectrum as a func- 
tion of frustration, E(f), The critical temperature gives 
a measure of the ground state energy, and in this respect 
we, as with Behrooz et al., find that the deep energy 
minima occur when it is simultaneously true that the 
magnetic flux through the large and small plaquettes is 
approximately an integer multiple of the flux quantum. 
We show that a model that treats the plaquettes indepen- 
dently describes these energy minima quite well. Consid- 
ering the set of energy minima in £"(/), we find that a 
quasiperiodic function emerges, with this quasiperiodic 
property being dictated by the incommensurate areas 
and not the aperiodicity. 

The structure of the paper is as follows. In Sec. II we 
introduce the Independent Plaquette Model. We begin 
our study of the aperiodic ladders in Sec. III. Conclusion 
comes as Sec. IV. An Appendix complements this study 
in that it provides a model with cosine interaction, having 
an analytic fractal solution for E(f). 

II. THE INDEPENDENT PLAQUETTE MODEL 
(IPM) 

Consider first a Josephson junction ladder in a uniform 
perpendicular magnetic field. The hamiltonian for this 
system is 

H = - Ej^[cos(9 lyk - 0i,fc+i) + cos(9 2 ,k ~ 2 ,k+i) + 

k 

+ cos(6i, k - 9 % h - 2nfx k )], (2.1) 

where 9\, k and 02, k are the phases of the superconducting 
order parameter along the upper and lower sides of the 
ladder, respectively. We have chosen the gauge of the 
vector potential so that it has zero component along the 
horizontal junctions. Then the effects of the applied field 
are represented by Ay = f? A(r)-dl — 2irfxk where x k 
is the x-coordinate of the /cth vertical junction (a rung of 
the ladder). The frustration factor / is the flux through 
a plaquette of unit area, measured in units of the flux 
quantum. 



FIG. 1: Two kinds of plaquettes with different areas that are 
part of a Josephson junction ladder. The arrangement of such 
plaquettes in the ladder can follow an aperiodic sequence. 



value (corresponding to the peaks in T c ) for values of / 
that can be understood by imposing the condition that no 
plaquette be frustrated^ We will show that the "Inde- 
pendent Plaquette Model" is also useful for the aperiodic 
situation. 

First, let us consider a single isolated plaquette con- 
taining four junctions; furthermore, let us choose the 
gauge so that A+j has the same value for each successive 
junction as we go around the plaquette. Then we will 
have A = irfS/2, where S is the area of the plaquette, 
and the energy becomes 

H = - cos{6 x —02— A) — cos{8 2 — 83 — A) 

- cos(9 3 -6i-A)- cos{9 i -61 — A), (2.2) 

where 9i are the phases at successive sites of the ring. 
Analysis shows that the minima of H occur for 9i such 
that 9i~9i + \ — nn/2 where n is an integer; the minimum 
energy per junction is given by 

E = -cos(^(fS-[fS})), (2.3) 

where [x], throughout this paper, denotes the integer 
closest to x. It is readily seen that the minimum en- 
ergy per junction occurs when fS is an integer, and that 
it is periodic in / with period 1/S. Near a minimum, 
the cosine function is relatively insensitive to variations 
in the 9i away from the optimum values. This suggests 
that even in an array, when one plaquette is close to its 
energy minimum it will not be greatly perturbed by its 
links to other sites, nor will its presence greatly perturb 
the other rings. Then we may attempt to approximate 
the ground state energy per junction for an interacting 
system by averaging the plaquette energies Eq. (|2.3[) . 

In the present case, the ladder has two types of plaque- 
ttes; L with area Al, and S with area As- Let N Ladder be 
the number of plaquettes in the ladder, Nl the number 
of plaquettes of type L, and Ng the number of plaquettes 
of type S. Then within the IPM the energy per junction 
of the ladder becomes 

Eipm = {N S E S + N L E L ), (2.4) 

M Ladder 

where Es{El) is the energy per junction for a plaquette 
of type S(L). From 12.31 we can write: 

E s = -co S (~(A s f-[A s f})), (2-5) 
E L = -cos(^(A L f-[A L f})). (2.6) 

It will prove convenient to choose As to have unit area 
(which sets the scale for /), and define a = Al/As and 
v = Nl/Ns- Then the average spacing between vertical 
junctions (the rungs of the ladder) is 



For a periodic ladder, the ground state energy is a peri- 
odic function of /. The energy takes on the unfrustrated 



a = 



1 + vo. 



(2.7) 
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and we can write Eqs.(2.5 - 2.7) as 

Eipm = t^—( E s + vE L ), (2.8) 
1 + v 

E s = -cos{\{f-[f])), (2-9) 

E L = -cosfyaf - [af])). (2.10) 

We will not claim that the IPM is anything more than 
a crude approximation. Neighboring plaquettes are al- 
lowed to assign different currents to the rung they share; 
not only is this ignored, the rung energy is counted twice 
in the averaging! However, it is easy to use, and we will 
see that its predictions are not too different from the nu- 
merical simulations (to which it seems to serve as a lower 
bound). We will be particularly interested in the values 
for / for which the energy per junction is close to the 
absolute minimum —Ej, and for these cases the super- 
current on every junction is small, and the predictions of 
the IPM become accurate. 

III. STUDY OF APERIODIC LADDERS 

Here we compare a numerical study of the spectrum of 
the ground states of various frustrated Josephson junc- 
tion ladders with the predictions of the IPM. Each ladder 
is a sequence of two kinds of plaquettes, L and S. This 
enables us to distinguish the importance of two proper- 
ties of the ladder: the ratio of areas of the two kinds of 
plaquettes, and the order or particular arrangment of the 
plaquettes in the ladder. 

The ground state energy E(f) is found numerically us- 
ing the Monte Carlo (MC) simulated annealing routined 
A candidate ground state configuration is what is ob- 
tained at the lowest temperature attempted, which we 
have taken to be T — 0.001. (By this we mean kpT/ Ej = 
0.001.) 

We can detect any kind of order that might exist in 
ground state energy as a function of frustration, by cal- 
culating the "discrete power spectrum"— for E(f) known 
for N values of / in a range —F < f < F (E(f) is an 
even function of /, so the range can be doubled). The 
Fourier transform of E(f) is 

E{u) = ^ J E(f)exp(i2nfcj) 
1 W_1 

w — E(nF/N)exp(i2nk/N) (3.1) 

71=0 

where u) = k/F with k = 1 - N/2, N/2 - 1. Then 
the discrete power spectrum is 

S u = (\E bJ \ 2 + \E^\ 2 ), (3.2) 

To the extent that E(f) can be represented as a sum of 
periodic components, S(uj) will have peaks at the corre- 
sponding frequencies. 



For better comparison, we define an incomplete sum 
over frequencies as 

3 

IS u= j/p = ^ S k . (3.3) 
fc=i 

When at a frequency j the power spectrum has a peak, 
the incomplete sum of the discrete power spectrum has a 
jump. Therefore, jumps in IS U indicate peaks in S^, with 
the size of the jump indicating the summed weight of the 
peak. In the results to be presented we used an increment 
of 1/13 in frustration and 256 data points, except for the 
case of crFibonacci where we had an increment of 1/55 
and 5500 points. 

The power spectrum was calculated in MATLAB by 
the FFT method^ 



A. r Fibonacci Ladder 

We start with a ladder made of two types of plaque- 
ttes, with area ratio r and arranged according to the Fi- 
bonacci sequence; we call this the rFibonacci ladder. The 
plaquette areas are A L = t = (1 + v / 5)/2 and A s = 1. 
Denoting the structure in its nth step of construction by 
U n , the structure is constructed recursively by the rule 
U n +2 = U n +\ + U n , where summation means adjacent 
placement: to get the structure in step (n + 2), put the 
structure in step (n + 1) first, and that of step n, to its 
right. This generates the sequence 

S, L, LS, LSL, LSLLS, LSLLSLSL, ... . (3.4) 

In the limit of an infinite structure this defines a qua- 
sicrystalline ladder with v = t. 

The quasiperiodic grid defining the sequence of long 
and short spacings for our rFibonacci ladder is given byAS 

x n =n + 5 + r- 1 [T- 1 n + 13}. (3.5) 

where S and /3 are arbitrary constants, chosen to be 5 = 0, 

Within the IPM each plaquette minimizes its energy by 
itself; in this way the expression for the energy becomes 

Eipm = --^co S ^(f-[f}))--^cos(^(rf-[rf})). 

1 + T Z 1 + T i 

(3-6) 

Eq. (|3.6I) is composed of two terms; the first has period 
A/ = 1, while, the second has period A/ = 1/r. Hence 
this function attains its minima for those / which are 
close to integers and which bring rf close to an integer. 
These include the Fibonacci numbers 2, 3, 5, 8, 13, 21, but 
there are other values of / (/ = 6.11,6.88, 11.08, ... that 
also meet this criterion. Even though the energy as a 
function of frustration is not periodic, it shows order in 
its structure. For example, the energy spectrum between 
the special values of / that minimize frustration is nearly 
symmetric about the midpoint; e.g. E(5 + x) ~ E(8 — x). 
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The numerical calculation of the ground state energy 
is given in Fig. 2 along with the predictions of Eq. (|3.6p . 
It is not surprising that Ejpm < EmC: because the IPM 
neglects the interaction between the plaquettes. How- 
ever, most features of the numerical curve are reproduced 
by the IPM. We can quantify this comparison by calcu- 
lating the incomplete integral of power spectra, shown in 
Fig. 3. 
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FIG. 2: Energy for the rFibonacci ladder. The area ratio is 
r = (1 + v / 5)/2 and the order of the plaquettes is determined 
by the Fibonacci sequence. The upper curve is the ground 
state energy as determined by simulated annealing; the lower 
curve is the result of the Independent Plaquette Model (IPM) 
approximation. 
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FIG. 3: Incomplete sum of the Fourier spectrum of E(f) for 
the rFibonacci ladder using the numerical minimization (top) 
and the IPM (bottom). The top figure shows jumps at ui — 
t — 1, 1,t, t + l.This shows that the energy is quasiperiodic 
as a function of /, with periodicities A/ = t~ 2 ,t~ , 1, and 
r, and weaker periodicities at various integer combinations of 
unity and r. Figures are shifted relative to each other, for 
clarity. 



Both E(f) contain significant periodicity correspond- 
ing to the frequencies uj = 1 and w = r. While the IPM 
energy expression is a sum of two periodic functions with 
incommensurate periods, the numerical result indicates 
that this function is quasiperiodic too, as it also includes 
the sum and difference of the frequencies in the power 
spectrum: r + 1 and r — 1. This is also what is found 
for the tight binding hamiltonian model on a Penrose 
latticed 



B. crFibonacci ladder 

We have seen that the quasiperiodic nature of the se- 
quence and the special ratio of the areas both play roles 
in the power spectrum. To distinguish the role of the 
area ratio from the effects of the sequence we will keep 
the Fibonacci order for the arrangement of the S and L 
plaquettes, but choose the areas to be in the silver ratio: 
Ai, = (7 = 1 + y2 and As = 1. The general formula for 
generating this ladder is now-i£ 

x n = n + S+[-+(3]V2, (3.7) 

T 

with 8 and (3 being arbitrary, chosen here to be 6 = 0, 
P =-1/2. 

For this ladder the energy within the IPM is expressed 
as, 

Eipm = -^C0s£{f-[f]))--?—cos£{af-[af])). 

(3.8) 

Figure 4 shows the result of a simulated annealing 
Monte Carlo study of the ground state energy for a lad- 
der with 144 plaquettes, along with the predictions of 
the IPM, Eq. (3.5). The properties mentioned for the 
rFibonacci ladder hold here too: there are minima near 
/ = 2, 3, 5, 7, 10, and 12 because for these values erf is 
also close to an integer. The comparison of the Fourier 
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FIG. 4: Energy for the oFibonacci ladder. The area ratio is 
a = 1 + V% and the order of the plaquettes is determined by 
the Fibonacci sequence. The upper curve is the ground state 
energy as determined by simulated annealing; the lower curve 
is the result of the Independent Plaquette Model approxima- 
tion. 

spectra for the two energy plots, Fig. 5, indicates that 
the dependence of the energy on / is dictated by a and 
not by the specific way the plaquettes are arranged. The 
ratio of the areas causes the frustration, and the energy 
minima occur for the values of / that give nearly integer 
flux through each plaquette. The quasiperiodic nature 
of E{f) indicates that there is a fine structure to it that 
corresponds to a 'many tile effect,' pointing to the pres- 
ence of vortex lattices commensurate with the underlying 
structure 4^ The IPM cannot predict this effect, and in 
Fig. (5) we have peaks that are not reproduced by the in- 
dependent plaquettes approximation, just as in the case 
of Fig. (3). 
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FIG. 5: Similar to Fig. 3 but when E(f) is for the aFibonacci 
ladder. The jumps at uj — a — 1, 1, a, a + 1 show that the 
energy is quasiperiodic as a function of /, with periodicities 
A/ = (<t+1) -1 , c -1 , 1, and (<j— and weaker periodicities 
at various integer combinations of unity and a. Figures are 
shifted relative to each other for clarity. 



C. Other r Ladders 



FIG. 7: Periodicity of E(f) for the rRandom ladder, as rep- 
resented by the incomplete integral of the power spectrum 
(similar to Fig. 3). It shows jumps at u = t — 1, 1, r, r + 1. 
This shows that the energy is quasiperiodic as a function of / 
with periodicities A/ = r~ 2 ,r _1 , 1, and r, and weaker peri- 
odicities at various integer combinations of unity and r. The 
characteristics in this figure are nearly identical with those of 
the other two rLadders. Figures are shifted relative to each 
other, for clarity. 



To further study the importance of plaquette order, 
we have considered three ladders with v — 1 and a = r, 
having different arrangements of the S and L plaquettes: 
random order, the Thue-Morse sequence^ and a peri- 
odic lattice. 

For the random ladder, the results for a simulation 
using 144 plaquettes is shown in Fig. [5] According to 
our simulated annealing study, there is less fine struc- 
ture in the random ladder than predicted by the IPM . It 
remains true that the IPM consistently gives a slightly 
lower energy that the simulated annealing, and that it lo- 
cates the energy minima. The incomplete integral of the 
power spectrum of the rRandom E(f) is shown in Fig. 
[7J the other two rLadders mentioned here give similar 
results. 




FIG. 6: Energy for the rRandom ladder. The area ratio is 
r = (1 + v / 5)/2 and the order of the plaquettes is determined 
by the random sequence. The upper curve is the ground state 
energy as determined by simulated annealing; the lower curve 
is the result of the Independent Plaquette Model (IPM) ap- 
proximation. 

We have also studied the rPeriodic lattice (a periodic 
sequence of SL, with areas having ratio r). as well 
as one having the Thue-Morse order (built recursively 
by appending the complement of a sequence to itself: 
SL, SLLS, SLLSLSSL, .... ). Our results are summa- 
rized in Table U for the numerical results and Table |TT] 
for the IPM results. We see that all lattices have peaks 



at u> = 1, a and a ± 1, and this shows that E(f) is a 
quasiperiodic function with main peaks at u> = 1 and a. 
From these results we separate out the role of the areas of 
plaquettes from that of the order in lattice. rFibonacci 
and crFibonacci lattices have the same order but with 
different areas of plaquettes. The results for these two 
lattices show that the position of the main peaks is de- 
termined by the areas of plaquettes. 

We can further justify the results by a quantitative ar- 
gument that starts with the IPM model. From the IPM 
model we know that the main peaks in S(u>) are propor- 
tional to the square of the density of S and L. These are 
just the coefficients of E$ and El, in the expression for 
energy (Eqs. 2.8-2.10). Then the peak heights depend 
only on v; this is why the last three rows of Table II (for 
different v = 1 models) have identical entries for the cor- 
responding peaks. In the same way, we can explain the 
added peaks in the MC data of Table I in terms of the 
correlations between plaquette types. We define £>ls to 
be the probability that the two types of plaquettes ap- 
pearing adjacent to each other: the ratio of all possible 
LS and SL pairs to the total possible number of pairs in 
the array. This is correlated with the height of the peaks 
at oj = a + 1, as shown in Table III. For the r ladders, the 
peak height is proportional to D^s to good approxima- 
tion. There is a similar explanation^ for the amplitudes 
of the peaks at u) = r — 1 that uses the LS combination as 
its basic unit, so that two-plaquette interactions account 
for this case, too. 

In summary, our study of the one dimensional lat- 
tices shows that E(f) is a quasiperiodic function with 
the main frequencies at the areas of the plaquettes. This 
property is an extension of that of the periodic behavior 
for the square latticed We can also observe an approxi- 
mate mirror symmetry in the dependence of E on /; by 
that we mean, there is a set of values /„ for which the 
energy is very close to its unfrustrated value, and be- 
tween these values E(f) has a mirror symmetry about 
/ = (fn + /n+i)/2. This can again be viewed as an ex- 
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rFibonacci 


t-1 (0.653) 


1 (2.360) 


r (12.14) 


r+1 (3.518) 


crFibonacci 


1 (3.457) 


cr-1 (0.397) 


a (14.0) 


cr+1 (3.249) 


rperiodic 


t-1 (0.565) 


1 (7.132) 


T (8.111) 


r+1 (3.899) 


rThue-Morse 


t-1 (0.422) 


1 (10.16) 


r (9.387) 


r+1 (2.334) 


rRandom 


t-1 (0.5285) 


1 (7.138) 


r (8.164) 


r+1 (1.455) 



TABLE I: Peak locations and heights for different lattices, 
according to the simulation. For each lattice four main peaks 
given; number in parentheses gives height of peak (we multi- 
ply it by 1000 for better representation). 



rFibonacci 


1 (4.208) 


r (11.008) 


crFibonacci 


1 (4.208) 


a (11.008) 


rperiodic 


1 (7.2) 


t (7.2) 


rThue-Morse 


1 (7.2) 


t (7.2) 


rRandom 


1 (7.2) 


t (7.2) 



TABLE II: Peak locations ui and heights (in parentheses) for 
different lattices, according to the IPM. For each lattice two 
main peaks given; number in parentheses gives height of peak 
(we multiply it by 1000 for better representation). 



tension of the mirror symmetry for the square lattice. 



IV. CONCLUSION 

The independent plaquette model treats the plaque- 
ttes as if they were independent of each other, and at 
the deep energy minima, this nearly coincides with the 
energy minimum when such an approximation is lifted. 

By using the numerical minimization of Hamiltonian 
12.11 we find some interesting results about the structures 
of ground state energy E(f). We see that E(f) has two 
main frequencies and their ratio equals the ratio of areas 
of the two plaquettes. There are other frequencies in the 
spectrum of E(f) equal to sum and difference of the two 
main frequencies, indicating that E(f) is a quasiperiodic 
function and not just a sum of two periodic functions with 
incommensurate periods. The order of the lattice affects 
the amplitudes of different frequencies in the spectrum of 
E(f). Study of the energy power spectrum gives a role to 
the lattice ordering because this affects the correlations 



rFibonacci 


2/r 2 


6.028 


crFibonacci 


2/r 2 


5.567 


rperiodic 


1 


5.899 


rThue-Morse 


2/3 


5.252 


rRandom 


1/2 


5.820 



TABLE III: The ratio of all possible LS and SL pairs to 
the total possible number of pairs in the array is defined as 
Dls, and given in first column. The second column gives the 
ratio of S(a + 1) to D\ s (we multiply it by 1000 for better 
representation) . 



between S and L, which goes beyond the IPM. However, 
this is only a short range effect; we have no evidence that 
the long range order in the array plays any role. 

If we look at the behavior of T c (f), which roughly goes 
hand in hand with the ground state energy, the linearized 
mean-field theory provides an analytical approach to the 
problem^ The results of this model show that T c (/) is a 
quasiperiodic function^ In their respective Fourier spec- 
tra, we observe peaks on the same special combinations 
of the main frequencies, and this is an indication of the 
particular order present in the lattice. 

The independent plaquette model provides an excellent 
approximation for E(f), following nearly all its features, 
while apparently giving a lower bound for it. 
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APPENDIX: IPM and the J model 

In a study with similar motivation, Grest, Chaikin and 
Levine 2 ^ (GCL) numerically studied the ground states 
of a one dimensional quasiperiodic array with inflation 
symmetry. Although their model was inspired by the 
Josephson network, it differed from it in two ways: 

*They simplified the form of the interaction, in effect 
replacing cos(x) by 1 — i(x — 2im) 2 (where n is the inte- 
ger that gives the minimum value for this expression 2 !) . 
In practice the arguments of the cosine function are small 
(mod 2tt) , so that this is a good approximation; it some- 
what changes the energy cost of introducing a vortex, 
but this is unimportant since the number of vortices is 
set by the applied field. They refer to this as the "current 
squared" (J 2 ) model. 

*They assumed that no current is carried along the pe- 
riodic wires* 2 ^ This forces the phases along these wires to 
be all the same (in a certain gauge), and implies that all 
plaquettes along a column either all contain a vortex or 
all do not. This reduces the mathematics problem to one 
dimension, so that Griffiths and Florfa 2 ^ (GF) were able 
to solve the model exactly. In their numerical study, GCL 
found that there are deep minima in E(f) which they in- 
terpreted as due to a kind of commensuration between 
the current pattern and the underlying lattice. Similar 
results were found in experimental studies by Chaikin et 
ali^Si For the J 2 model, GF showed that E(f) takes on 
the value 1 — 7r 2 /3 (the averaged value of the approxi- 
mated cosine function) for almost all /, but that there 
is a countable set of irrational values of / at which E{f) 
takes on a lower value. 

We have some misgivings about the choice of model. 
The way that GCL reduce it to one dimension does some 
violence to the physical situation. A column of vortices 
is unstable, because the vortices repel each other. We 
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would expect that the ground state configurations would 
have the vortices more uniformly spaced in two dimen- 
sions. Having no current on the periodic wires also means 
that the current is the same on all of the nonperiodic 
wires and does not vary along it. We will show that with 
this very strong constraint, it is not necessary to make 
the J 2 approximation: the model can be exactly solved 
for the cosine interaction (with very little difference in 
the results). 

The GCL model is defined by the energy function 
H = -EjJ~] [cosjOjj - )+cos(6 i> j - 9 iij+ i - 2nfxj )] . 

(A-l) 

As an approximation, GCL assumed 9ij to be indepen- 
dent of i, which implies that there is no current along 
the "horizontal" links corresponding to the index i. This 
makes the first term a constant which will not be included 
in what follows. Minimizing the energy with respect to 
9j gives 

sin(0j — Oj + i — 2nfxi) — sin(9j-\ — 8j — 2Trfxi), (A-2) 

which implies 6j — 9j + i = 9j-\ — 9j (mod(27r)). Then 
there is just one undetermined variable: the difference 
in phase b between adjacent columns. Specification of 
b determines all of the currents and link energies in the 
form 

1 N 

E=- lim — Ej Vcos(6- 2tt/x,-)- (A-3) 

j=l 

This is not the same as the model that we have been 
discussing, but the results can be compared to those of 
the IPM. Applying the IPM to the two-dimensional ar- 
ray again leads to Eqs. (2.8 - 2.10) since it decouples 
neighboring plaquettes. 

Eq. (|A-3[) can be evaluated exactly following the meth- 
ods of GF. They define r{y) to be the distribution of the 



values yj = fxj mod(l), so that the energy can be cal- 
culated from 

E = -Ej [ cos(2ny)r{y - b)dy. (A-4) 
Jo 

They consider sequences {xj} that can be written in the 
form Xi = a[i+g(iuj)] where g{x) is the periodic extension 
(with period 1) of the function g(x) = 72;. These include 
the periodic, r, and a ladders — for example, the r ladder 
corresponds to a — 3 — r, 7 = l/^/(5), and w = r — 1. 
They show that for almost all values of /, r{y) is the 
uniform distribution r(y) — 1 (so that E = 0), but that 
r(y) has nontrivial form for af = (I + muj)/n, where I, 
to, and n are integers and a sets the scale. The effects 
for n > 1 are relatively small and will not be considered 
further. Then 



r(y) = (K+1)/k for - s/2 < y < s/2, (A-5) 

= K/k otherwise. (A-6) 

Here, we have s — n — K , with k = |(1 + 7cj)m +7?|, 
and K the largest integer less than k. The corresponding 
energy is2£ 

E = -Ej\sin(iTK)/TTK\. (A-7) 

Because GF follow GCL making the approximation 
cos(x) = 1 — x 2 /2, they obtain a slightly different result. 

It is readily seen that the most negative values for 
Eq. (|A-7p occur when k is small, i.e. when m/l ~ 
|1 + 7cj|/|7|. The dependence of the energy on / de- 
scribed by Eq. (|A-7[) is a discontinuous fractal function — 
a forest of vertical lines, rather than the smooth curves 
we find numerically for the Josephson ladder. We believe 
this to be a real difference between the two models. In 
any case, the two results agree on the positions of the 
deepest minima. 
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